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We discuss the slow relaxation phenomenon in glassy systems by means of replicas by constructing 
a static field theory approach to the problem. At the mean field level we study how criticality in 
the four point correlation functions arises because of the presence of soft modes and we derive an 
effective replica field theory for these critical fluctuations. By using this at the Gaussian level we 
obtain many physical quantities: the correlation length, the exponent parameter that controls the 
Mode-Coupling dynamical exponents for the two-point correlation functions, and the prefactor of the 
critical part of the four point correlation functions. Moreover we perform a one-loop computation in 
order to identify the region in which the mean field Gaussian approximation is valid. The result is 
a Ginzburg criterion for the glass transition. We define and compute in this way a proper Ginzburg 
number. Finally, we present numerical values of all these quantities obtained from the Hypernetted 
Chain approximation for the replicated liquid theory. 



Contents 



I. Introduction 2 

II. An illustration of the main results of this paper in the simpler case of a standard 

ferromagnetic transition 3 

A. Landau theory 3 

B. Gradient expansion 4 

C. From the microscopic Hamiltonian to a field theoretical formulation 5 

D. Ginzburg criterion 6 

III. Dynamical heterogeneities and replicas: definitions 7 

A. Two-point functions: the dynamical order parameter 8 

B. Four point functions: the correlations of the order parameter 9 

C. Gonneetion between replicas and dynamics 9 

D. The replicated free energy 11 

IV. Landau expansion of the free energy around the glassy solution 13 

A. Free energy for a uniform field 13 

B. Gomputation of A 15 

C. The mass matrix and the zero mode 16 

D. Analysis of the cubic terms and computation of A 18 

V. Gradient expansion 19 

A. Fourier transform 19 

B. Spectrum of the mass matrix 21 

VI. Ginzburg criterion 23 

A. Projection on the zero mode 23 

B. ip^ theory in random field 23 

VII. Explicit calculations in the Hypernetted Chain (HNC) approximation 26 

A. Derivatives of the free energy in HNG 26 

B. The HNG mass matrix 27 

C. Expression of A in HNC 28 



2 



D. Computation of /i and cr in HNC 28 

E. Summary of the numerical results 29 

VIII. Conclusions 31 

Acknowledgments 32 

A. Double Counting Problem on the Lp^ theory 32 

References 33 



I. INTRODUCTION 

Much progress in the recent understanding of glassy relaxation of supercooled liquids has come from the study of 
dynamical heterogeneities [1]. Space-time fluctuations of the density field result in a distribution of regions of different 
mobility with typical size that grows upon decreasing the temperature, and that persist for time scales of the order of 
the relaxation time. Remarkably, the present theories of glassy relaxation are able to predict the qualitative features 
of this effect. Many progresses have been achieved neglecting activated processes. Both theories based on dynamics 
(Mode Coupling Theory [2]) and on constrained equilibrium (e.g. the molecular liquid theory [3]) agree in a critical 
instability of the liquid phase at a finite transition temperature. As far as universal aspects are concerned, a first 
theoretical insight originally came from spin glass theory, that suggested to look at the critical behavior of four point 
density correlation functions [4, 5]. In the context of p-spin models where both the equilibrium and the dynamic 
approach can be pursued exactly, it was possible to predict the qualitative features of dynamical correlations growth 
and associated dynamic criticality observed in numerical simulations [6] to well defined features of static correlations 
in constrained equilibrium. Beyond these schematic models a direct connection between equilibrium and dynamics is 
more difficult to make. However the static hallmarks of dynamical criticality are generically present whenever one has 
a one step replica symmetry breaking transition (or mean-field random first order transition [7, 8]), as it is found e.g. 
for liquids in the Hyperncttcd Chain (HNC) approximation [9]. At the level of dynamical liquid theory, the growth 
of correlations has been found with generality within Mode Coupling Theory [10] through a one-loop diagrammatic 
expansion. Despite these important results, in order to progress further it is necessary to build a field theoretical 
description of dynamical fluctuations capable in principle to go systematically beyond the zero-loop Gaussian level. 
Unfortunately, the dynamic approach is rather problematic in this respect for several reasons. The MCT is not a 
self-consistent theory as it requires the input of the static structure factor. Moreover, the equation for the dynamic 
correlator does not follow from a variational principle. Clearly this is not the ideal starting point to build up a 
perturbative computation. In addition an expansion around mean field in dynamics is prohibitively difficult even in 
the simplest cases. 

A static formulation would be therefore very useful to make additional progresses. The connection between statics 
and dynamics is based on the idea that the emergence of slow dynamics is due to the appearance of long lived 
metastable states [11, 12]. Within Mode-Coupling Theory and similar dynamical theories of glass formation, one 
can distinguish two separate dynamical regimes: the so-called /3 regime, that corresponds to the long time dynamics 
inside a metastable state, and the a regime, that corresponds to transitions between different metastable states. A 
general theory of fiuctuations in the /3 regime based on the replica method was proposed in [13] on the basis of general 
symmetry considerations. It was found that the dynamical transition is in the universality class of a cubic random 
field Ising model. Moreover in [14] it has been shown at the level of schematic models that the same theory can be 
used to evaluate the Mode-Coupling dynamical exponents and that they are related to the amplitude ratio between 
correlation functions that are in principle measurable. At the Gaussian level, we can consider this phenomenological 
theory as the static Landau theory of the glass transition (see [15] for a discussion of a dynamical Landau theory 
of the glass transition). In summary, the advantage of the approach discussed in [13] is that one can compute the 
long time behavior of the dynamical correlations in the (3 regime starting from a completely static replicated (and 
constrained) Boltzmann measure [16, 17]. The first consequence of this is that all the calculations simplify drastically. 
Moreover one can look at the dynamical transition through a static measure and one can see systematically where 
the mean field prediction is valid and where it fails. 

The aim of this paper is to construct a framework to study the glass transition in the /? regime by following as 
closely as possible the standard treatment of critical phenomena using field theoretical tools. Furthermore, we want to 
obtain the replica field theory of critical fluctuations starting directly from the microscopic grand-canonical expression 
for the replicated partition function for the liquid, in such a way that the couplings appearing in the field theory can 
be computed starting from the microscopic potential. This program can be achieved through an analysis of the soft 
modes that appear at the dynamical transition and that are responsible for the criticality. By computing these soft 
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modes we can focus on fluctuations that are along them: in this way we construct a gradient expansion for the field 
theory of the critical fluctuations and we compute in full details the Gaussian correlation functions, thus obtaining 
the critical part of the long time limit of the dynamical four point functions in the /3 regime. Then we look at the 
corrections to the Gaussian theory and we introduce a Ginzburg criterion in order to see where the Gaussian theory 
is valid. The Ginzburg criterion can be used in two ways. On the one hand it gives the upper critical dimension for 
the model, on the other hand it provides a measure of how much one has to be close to the transition line in order to 
see the non Gaussian fluctuations that cannot be treated by the mean field approach. 

A short account of our findings appeared in [18]. This paper is organized as follows. In Section II we review all 
the line of reasoning we just discussed for the standard Ising model. Then, in Section III we discuss how to obtain a 
replica description for the dynamical correlation functions in the j3 regime. In this way we rephrase the problem from a 
dynamical one to a standard static computation. In Section IV we discuss the expansion of the free energy around the 
glassy solution. By studying the Hessian matrix we identify the soft modes that appear at the dynamical transition 
and we compute the expression for the exponent parameter A that can be related to the dynamical Mode-Coupling 
exponents a and h that describe the critical slowing down of the two point correlation functions at the dynamical 
transition. Then, in Section V we show how we can perform a gradient expansion of the replica field theory in order 
to study the long distance physics and in Section VI we use this effective theory to compute where the Gaussian level 
computation fails by introducing a Ginzburg criterion for the dynamical transition. Up to that point we will remain 
completely general. The only assumption that we will make is that the replicated system has a glassy phenomenology, 
namely that the replica structure of the two point density function is non trivial below a certain dynamical transition 
point. In Section VII we report concrete calculations in the framework of the replicated HNC approximation by giving 
explicit expressions for all the couplings and masses of the effective replica field theory starting from the microscopic 
potential and we show the numerical results for several physical systems. 

II. AN ILLUSTRATION OF THE MAIN RESULTS OF THIS PAPER IN THE SIMPLER CASE OF A 

STANDARD FERROMAGNETIC TRANSITION 

The aim of this paper is to build a theory of the glass transition by following closely the first steps of the standard 
field theory formulation of critical phenomena. Namely, we want to start from a Landau theory (including microscopic 
parameters) and deduce from it a set of mean field critical exponents. This is done first by studying the behavior of the 
uniform order parameter in the mean field theory, and then considering a gradient expansion for slowly varying order 
parameter to compute the correlation length. Finally, taking into account the interaction terms lead to a pcrturbative 
loop expansion that allow to establish the region of validity of mean field theory (hence the upper critical dimension) 
through a Ginzburg criterion. The aim of this paper is to repeat all of these steps in the more complex case of a 
system undergoing a glass transition. The final step would be of course to set up an epsilon expansion around the 
upper critical dimension using renormalization group methods. We will not make any attempt in this direction in 
this paper. For pedagogical reasons, we find useful to briefly describe how these steps are carried out in a simple 
ferromagnetic system before turning to the glass case. The reader should keep in mind that this section is just a short 
reminder of the main steps, for a reader already accustomed with the modern theory of critical phenomena. 



Suppose we consider a microscopic system undergoing a ferromagnetic transition. In the following we will consider 
the ferromagnetic Ising model with nearest neighbor interactions on a D-dimcnsional cubic lattice and a properly 
scaled coupling constant: 



Starting from the microscopic Hamiltonian, we can construct the free energy as a functional of the order parameter, 
the magnetization = ((Ti), as follows. The free energy (here we will ignore some factors of temperature by rcscaling 
some variables) in presence of an external magnetic field is 



A. Landau theory 




(1) 




(2) 



{<T.=±1} 
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Taking a Legendre transform [19], wc define 

r[0] = ^/i*0,-w[/i*] , (3) 

i 

or in other words 

g-r[0] ^ g-/9^^W+E.'''K-0.) ^ (4) 

{<T,=±1} 

where h* is the solution of '^^^^^ = 0^. The function r[0] is the free energy of the system as a function of the 
magnetization field. 

In order to detect the phase transition, we want to investigate the small behavior of r[(/)]. Let us make a crucial 
assumption, that r[(/>] is an analytic function of (j) around (f) = 0. This assumption is plain wrong in finite dimensional 
systems at the critical point and below. However, let us for the moment forget about this problem and proceed with 
our discussion. We can consider a uniform magnetization profile (j) and expand T[(j)\ at small (j). From the symmetries 
of the problem, we know that 

r[</)] = F|i™2^2 + |04 + ...| , (5) 

which is the celebrated Landau free energy (here V is the volume of the system). A very practical way to compute 
the coefficients ttiq and g is to perform a systematic high temperature expansion^ [20]. For example, at the leading 
order, for the D-dimensional Ising model with coupling constant J = 1/{2D) wc obtain TOq = 1 — /?. Adding more 
terms leads to an expansion of mg in powers of /3. 

The equilibrium value of </), called (/), is obtained by minimizing T[(j)\. From the high temperature expansion we 
obtain that uiq vanishes linearly at a given temperature Tc, in such a way that TOq cx T/Tc — 1 = e. Note that for 
instance in D = 3 it is enough to consider the cubic term in the small /3 expansion to obtain a fairly accurate estimate 
of Tc- When mg becomes negative, the magnetization becomes non-zero with (j) ^ |to§|^/^ ^ e^/^ which gives one of 
the critical exponents. The other critical exponents are obtained in a similar way. 



B. Gradient expansion 



The subsequent step is to compute the correlation length. This is done by considering a gradient expansion for a 
slowly varying magnetization profile, again under the assumption that the expansion is regular at small 0. One can 
perform a continuum limit to simplify the notations: we denote by 1^3(2;) the continuum limit of the spin field it,;, while 
= (v(a;)) is the local average magnetization. The Landau free energy becomes at the quadratic order: 

^ dx(j){x){-\/^ +ml)(f){x) . (6) 



The correlation function of the magnetization is given by 

G{x -y)^ {Lp{x)ip{y)) 



d(j){x)d(j){y) 



Hence at this order the correlation function is^ 

Gq{v) = — 2 ' 



Go(x) = (^(.t)^(O)) 



(7) 



(8) 



which is often called the hare propagator. This expression shows that the correlation length is ^ = 1/mo ~ e ^^"^ and 
the magnetic susceptibility is x oc Go(p = 0) ~ e~^. 



^ For the Ising model (1), the true expansion parameter is actually /3J, i.e. the temperature in units of the coupling constant. Because 
the latter has to be chosen equal to J = 1/{2D) to obtain a good limit D — > 00, the expansion parameter is /3/(2D). In other words, 
the high temperature expansion is also a large dimension expansion around the D = 00 mean field limit. 

^ The calculation is performed by using l/(p^ + mg) = Jg°° dte~^^ +™o)* and changing variable to y = t/x^ . Then 

/■oo 2 1 

Go(x) oc x^-° / dye"*™"") ^ ^2-D_^(^^^) ^ 

Jo 

When X ^ 1/mOi a saddle point calculation shows that f{z) ~ z^^~^'>/'^e~'' , from which Eq. (8) follows. 
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C. From the microscopic Hamiltonian to a field theoretical formulation 



The above analysis relies on the assumption that T[(j)\ can be expanded as an analytic function around (j) ~ 0. 
Although this is certainly true at the mean field level (where the Landau theory provides the exact result), this is 
not the case in finite dimensional systems, because critical fluctuations induce a singular behavior of T[(j)\ at small (f). 
Hence we now want to assess the limits of validity of the Landau theory by studying the effect of critical fluctuations. 

The problem is that the definition of F given in Eq. (4) is not very convenient to perform a systematic expansion 
in the fluctuations around the mean field theory, although the computation could be done in principle. It would be 
much more convenient to write the effective action as a functional integral over a continuous spin field (p{x): 



with the following requirements: 

I. The mean field approximation should correspond to a saddle point evaluation of the above integral, in such a 
way that at the mean field level T[(j)] = S* [</>]. For consistency, S[(p] must therefore have a Landau form: 



in such a way that at the mean field level we recover Eqs. (5) and (6). In this way we can include fluctuations 
around mean field by performing a systematic loop expansion of the functional integral. 

2. The bare coefficients uiq and g entering in S[ip] must be reasonable approximations to the microscopic coeffi- 
cients as deduced from the Hamiltonian, in such a way that the mean field approximation is already a good 
approximation, and that loop corrections improve systematically over it. In this way we can guarantee that the 
criterion of validity of mean field theory has a quantiative meaning for the original microscopic Hamiltonian 



So wc want to give an appropriate definition of the continuum spin field 1^9(2;) and the corresponding action S'[(y9(a;)] 
in such a way that the requirements above arc satisfied. 

There arc many recipes for such a construction. Probably the best one is given by the non-perturbativc renormaliza- 
tion group construction [21], in which one defines a functional F£[(^(a::)] by integrating the small-scale spin fluctuations 
on length scales smaller than £, see e.g. [21, Eq. (28)]. If one chooses a "coarse-graining" length £0 that is quite bigger 
than the lattice spacing, but still quite small with respect to the correlation length (which diverges at the critical 
point), the function F£[(^(a;)] is an analytic function of (p at small (p, because the singularity is only developed at the 
critical point for ^ — 00 [21]. Then, Eq. (9) is basically exact with S replaced by F^^,. One can then expand F^^ 
at small ip and use this as the bare action S in Eq. (9). Although this strategy can be generalized to the physics of 
liquids [22, 23], calculations are quite involved so we need to consider something simpler. 

An alternative and very convenient prescription is the following. Let us call Ffc the truncation at a finite order (3'' of 
the high temperature expansion of F, as given in [20]. We know that Ffc[^] is an analytic function of cj) for any finite fc, 
hence Ffe[0] cannot be a good approximation of T[(j)\ at the critical point, because we know that F[(/)] is not analytic: 
in fact the high temperature expansion is divergent at the critical point. However, we can assume that Ffc[(y5] gives 
a good approximation for S[ip]. Note that our two requirements arc satisfied by the prescription that S[(p] = Tk[(p]. 
In fact, for the first requirement, at the saddle point level we obtain r[(/)] = S'[0] = F^ [(/)], and we already know that 
for fc = 1 this is the correct mean field result, while for /c > 1 we will obtain an "improved" mean field result. For 
the second requirement, we have already mentioned that the coefficients of Ffe[</)] give, for large enough k, a good 
estimate of the microscopic properties of the model (e.g. the critical temperature). Furthermore, we can argue that 
the high temperature expansion, at a given order k, is only sensitive to local physics up to a scale £{k) that grows 
with k. Hence, truncating the high temperature expansion at a finite order in k should be morally equivalent to 
perform an integration over the microscopic fluctuations on a scale smaller than i{k). We will see that this procedure 
is easily generalized to the case of liquids where the high temperature expansion is replaced by the low-density virial 
expansion. 

We will then use the prescription S[ip] — Ti;[(p], expand S[ip] in the form of Eq. (10), and use it in the functional 
integral Eq. (9) to compute F[(/)] in a loop expansion around mean field. Loop corrections give some non-singular con- 
tributions to F, which were already in part taken into account in the bare action S[ip] = Ffe[(p] because it was obtained 
from the high temperature expansion; hence we might have some "double counting" of non-singular contributions 
related to the short range physics. This double counting problem is discussed in more details in Appendix A. Still, 
our aim here is to find a Ginzburg criterion that identifies the region where these singular loop corrections are small, 
and the mean field approximations remains correct: we find that if the Ginzburg criterion is formulated in terms of 
physical quantities, then double countings are irrelevant. This is shown in next section II D and in Appendix A. 




(9) 




(10) 
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D. Ginzburg criterion 



We can use the above construction to perform a loop expansion in the couphng and check whether fluctuations are 
small such that they do not spoil the main assumptions we made above on the behavior of F at small 0, hence they 
do not change the critical behavior of the system. We will follow closely the derivation of [24] . Our bare action is 

1 f .^^ ..f^M V72 , _2^,.,^^ , 9 



S[^] = ^J dx ^(.x)(-V^ + m'^)ip{x) + ^ J dxip^ix) . (11) 

Here we will need to consider explicitly the presence of an ultraviolet cutoff (which will be of the order of the scale 
^o(fc) mentioned above). The bare propagator is 

Goip) = , (12) 

and the one loop correction to the propagator is [25] 

We can consider the inverse propagator and write 



G-\p) = Go(p)-i + 1 1 (^Go(g) , (14) 

which can be seen either as a Dyson resummation of the "tadpole" diagrams, or as an inversion of the perturbation 
expansion to obtain directly the second derivative of the Legendre transform of the generating functional. Physically, 
G~^{p = 0) is the "renormalized mass" or inverse magnetic susceptibility: 

m\ = G-\p = Q) = ml + '-j (3^)^^ ^2 + ^2 = -o + f / (s^-)^ ,2 + ^2^ ' (15) 

where the last equality of course holds at first order in g. The replacement of toq by is needed, because the 
perturbation theory must be done at fixed ruR, i.e. at fixed distance from the true critical point, otherwise corrections 
cannot be small [25]. The critical point is defined by the condition that mj:, = 0, or in other words the susceptibility 
is divergent, hence at the critical point 



2 = _ff / dg 1 
2 J i2n)D q2 

We see that the shift of the critical temperature is divergent in the ultraviolet (UV divergent) for D > 2: indeed, this 
is a non-universal quantity and depends on the details of the UV rcgularization. Now if we define the distance from 
the critical point as 

2^.9 dq 1 

we can write Eq. (15) as 

-<('-lfwF7l?^)^ '''' 

This is the crucial relation that relates the inverse susceptibility to the distance from the critical point. The Ginzburg 
criterion is obtained by imposing that the one loop corrections do not change the mean field behavior t = rnj^. 
We can now distinguish two cases: 

• For D < A, the correction is UV convergent and infrared (IR) divergent. In this case, we can safely send the 
cutoff to infinity because the renormalized theory exists. We obtain 



5 /"^ d^g 1 \ _ 2 _9_ D-2 /'°° 1 D-l 1 

2 J (27r)^g2(g2 + ^2^); 2 {2n)D x^ix^ + 1) 



7 



and the integral over x is finite. We clearly see that because D < 4, the second term will be dominant over the 
first close enough to the critical point. By imposing that the first term dominates, we obtain the criterion in 
the form 

1 » g m^-^Cu = ge-'^Co = Gi , (20) 

where we used that in the mean field region the correlation length ^ — X/mn. Hence the Ginzburg number 

Gi = gCo is a universal constant in this case. This shows that loop corrections will always be relevant close 

enough to the critical point and gives a precise value of the correlation length at which they will become relevant, 
l/(Gi)i/(4-i5). 

• Instead, for -D > 4, the correction is UV divergent and IR convergent. In this case the Ginzburg criterion is 
non-universal and strongly dependent on the details of the regularization. For a fixed cutoff A, the integral is 
finite at m|j. = and the mean field behavior is always correct: 

However one loop corrections provide a strong renormalization of the coefficient relating t to m|j. Imposing that 
these conditions are small we obtain 

9 d^q 1 



1 ^ i / ^ z (22) 

2 J {2t:)D q^iq-i +m\) ^ ' 

This provides a condition on ruR for a given UV cutoff A. When the condition is satisfied the mean field 
calculation is not only qualitatively, but also quantitatively correct. Note that the integral is upper-bounded by 
its value in uir = 0. Then, if 1 ^ gCA^^^, the Ginzburg criterion is always satisfied and one loop corrections 
are small even at the critical point. Instead, if gCA^~* is of order 1 of bigger, then we obtain a non-trivial 
condition on mfl and one loop corrections are large close enough to the critical point. 



III. DYNAMICAL HETEROGENEITIES AND REPLICAS: DEFINITIONS 

The aim of this paper is to repeat the steps outlined in Sec. II in the case of a glass transition. As we will see, the 
calculation is in this case complicated by several problems: 

1. We will need to introduce replicas to define a proper static order parameter. 

2. The order parameter is in general not a real number (e.g. the magnetization) but a function g{x — y) that 
encodes the replica-replica correlations. Hence we will have to introduce a smoothing function to define a scalar 
order parameter q ''-^ J f{x)g{x). We will then have to show that the choice of the function f{x) is irrelevant. 

3. The glass transition is discontinuous in the order parameter, which jumps to a finite value at the transition. 
Hence the transition is not an instability of the high temperature solution, but rather a spinodal point where 
the low temperature solution first appear. Because of that, we need to control the effective free energy at values 
of the order parameter that are very far away from the high temperature solution. Keeping only a few terms in 
the high temperature expansion is not enough, and we will have to resum an infinite number of terms to obtain 
a good starting point for the mean field theory. 

4. Because the glass transition is akin to a spinodal point, the resulting effective action is a cubic theory. Hence 
the theory is not really defined (spinodals do not exist in finite dimension). This will not be a problem for the 
mean field and loop calculations, but wc expect it to be a serious problem if one wants to go beyond mean field 
and construct a systematic epsilon expansion (which we do not attempt here). 

In this section we will better explain the first two points: we will give some important basic definitions on criticality 
at the glass transition (as encoded by the so-called dynamical heterogeneities) and wc will show how the problem can 
be tackled using replicas. In Sec. HI A we introduce the basic dynamical order parameter of the glass transition, and 
in Sec. HI B its correlation function. In Sec. HI C we show how both quantities can be written as static correlations in 
a replicated theory. In Sec. HI D we set up the general form of this replicated theory and give some useful definitions. 
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Throughout this paper we consider a system of N particles in a volume V interacting through a pairwise potential 
v{r) in a D dimensional space. The basic field is the local density at point x and time t: 

N 

p{x,t)=J2S{x-x,{t)) . (23) 

i=l 

We will consider a generic dynamics that can be either Newtonian or stochastic, e.g. of Langevin type. In both 
cases, we will consider equilibrium dynamics, that starts from an equilibrated configuration of the system. It will be 
convenient to separate the dynamical average in two contributions [13]. A dynamical history of the system will be 
specified by the initial configuration of the particles {xi(0)}, and by a dynamical noise. For Newtonian dynamics, this 
noise comes from the initial values of the velocities, extracted by a Maxwell distribution; for stochastic dynamics, it 
comes from the random forces that appear in the dynamical equations. We will denote by (•) the average over the 
dynamical noise for a fixed initial condition; and by E[»] the average over the initial condition. Hence, for instance, 
the equilibrium average of the density will be p = E[(/5(x, t))]. 



A. Two-point functions: the dynamical order parameter 

The dynamical glass transition is characterized by an (apparent) divergence of the relaxation time of density 
fluctuations, that become frozen in the glass phase. The transition can be conveniently characterized using correlation 
hmctions. Consider the density profiles at time zero and at time t, respectively given by p{x,0) and p{x,t). We can 
define a local similarity measure of these configurations as 



C(r,t)= f dxfix)p{r + ^,t)p{r-^,0) 



a;,(t)+x,(0)\ ^, (24) 



Sir- , M fix4t)-x,m 



Here f{x) is an arbitrary "smoothing" function of the density field with some short range a, which is normalized in 
such a way that / dxf{x) ~ 1. As an illustration, let us choose f{x) = d{a — \x\)/Vd{a), where Vd{a) is the volume of 
a sphere of radius a, and suppose that a is much smaller than the intcr-particlc distance and that t is short enough. 
In this situation, f(xi{t) — Xj(0)) vanishes unless i = j, and we get 

C{r,t)^Y.^('-^'^^^) fiMt)~MO))- (25) 

i ^ ^ 

Therefore, C(r,t) counts how many particles that are around point r have moved less then a in time t and is often 
called "mobility" field. Alternatively, Eq. (25) can be taken as the definition of a self two-point correlation function. 
Different choices of f{x) lead to other correlations that have been used in different studies. We will show later that 
the choice of the function f{x) is irrelevant as far as the critical properties are concerned. 
Let us call 

C{t)^V~' J drn{C{r,t))]-p' (26) 

the spatially and thermally averaged connected correlation function. Typically, on approaching the dynamical glass 
transition Td, C{t) displays a two-steps relaxation, with a fast relaxation" occurring on shorter times down to a 
"plateau", and a much slower "a-relaxation" from the plateau to zero [2]. Close to the plateau at C{t) = Cd, one 
has C(t) ^ Cd + At^"" in the /3-regime. The departure from the plateau (beginning of a-rclaxation) is described 
by C{t) ^ Cd — Bt^. One can define the a-relaxation time by C(tq,) = C(0)/e. It displays an apparent power-law 
divergence at the transition, ^ \T — Td\~^ . All these behaviors are predicted by MCT [2], which in particular 
relates all these exponents to a single parameter A through the formulae 

T{l-a? ^ n^ + bf ^ 

r(l-2a) r(l + 25) ' ^27) 
_ 1 1 
^ 2^ ^ 26 ' 

and gives a microscopic expression of A in terms of liquid correlation functions [2]. In low dimensions, a rapid crossover 
to a different regime dominated by activation is observed and the divergence at Td is avoided; however, the power-law 
regime is the more robust the higher the dimension [26, 27] or the longer the range of the interaction [28]. 



9 



B. Four point functions: the correlations of the order parameter 



It is now well established, both theoretically and experimentally, that the dynamical slowing is accompanied by 
growing heterogeneity of the local relaxation, in the sense that the local correlations C{r,t) display increasingly 
correlated fluctuations when is approached [1, 5, 29, 30]. This can be quantified by introducing the correlation 
function of C(r,t), i.e. a four-point dynamical correlation 



G4(r, t) = E[(C(r, t)C{0, t))] - E[(C(r, t))]E[(C(0, t))] 



(28) 



This function describes the total fluctuations of the two-point correlations, and it decays as Gi{r,t) ^ exp(— r/^(t)) 
with a "dynamical correlation length" that grows at the end of the /3-regime and has a maximum £, = £,{t ^ To) that 
also (apparently) diverges as a power-law when is approached. MCT [2] and its extensions [10, 31-35] give precise 
predictions for the critical exponents. 

For later convenience, we can also consider a modified four-point correlation: 



Gth (r, <) = E ((7(r, t)C{{), t)) - {C{r, t)) ((7(0, t)) 



(29) 



This function describes the isoconfigurational fluctuations of the two-point correlations, i.e. the fluctuations due to 
the noise of the dynamical process at fixed initial condition. It describes the in-state susceptibility: indeed, the initial 
condition selects a typical glass state, which is then explored by the dynamics. A final average over initial conditions 
is taken to ensure that the initial condition is a typical one. 

For each of these correlations, we can define the corresponding susceptibility 



XA{t) = j drG4(r,0 =E 
Xth{t) = / drGth{r,t)^-E 




drC{r,t) 



- I drC{r,t) 




drC{r,t) 



- I drC{r,t) 



(30) 



C. Connection between replicas and dynamics 



The dynamical glass transition can be also described, at the mean field level, in a static framework. This has the 
advantage that calculations are simplified so that the theory can be pushed much forward, in particular by constructing 
a reduced field theory and setting up a systematic loop expansion that allows to obtain detailed predictions for the 
upper critical dimension and the critical exponents [13]. Moreover, very accurate approximations for the static free 
energy of liquids have been constructed [36] , and one can make use of them to obtain quantitative predictions for the 
physical observables. 

In the mean field scenario, the dynamical transition of MCT is related to the emergence of a large number of 
metastable states in which the system remains trapped for an infinite time. At long times in the glass phase, the 
system is able to decorrelate within one metastable state. The /3 regime is identified with the dynamics "inside a 
metastable state", while the a regime is identified with "transitions between different states". Hence we can write 
(introducing two new averages): 

iCir, t^oo))^J d./(.) (p(. + |))Jp(.-|))^^^ , ^^^^ 

E[(C'(r, t ^ oo))] = J dxfix){p (r + |) )^ (p (r - |) )^ ■ 

In fact, if one performs a dynamical average at fixed initial condition, the system is trapped in a single metastable that 
can be explored, and at long times the dynamical average can be replaced by the average (•)m in the metastable state 
selected by the initial condition. The average over the initial condition then induces an average over the metastable 
states with equilibrium weights, that we denoted by an overline. 
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For the four-point functions wc obtain 

G4(r,t ^ J dxdyfix)f{y) [(/5(r - a;/2)p(-y/2))„(p(r + x/2)p(y/2)) J 

- E[(C'(r,t -> oo))]E[(C'(0,t cx)))] , 

f (32) 

Gthir, t oo) = / dxdy f{x)f{y) [(/)(r - x/2)p(-y/2)),„(p(r + a;/2)p(y/2)) J 



dxdy [{p{r - x/2)U{pi-ymMr + x/2)p(y/2))„] . 

The reason for the particular structure of the second term of Gth is that the densities at time come from the same 
initial condition and are therefore correlated, but they then evolve separately and therefore the two densities at time 
t are uncorrelated. 

The above structure suggests that the dynamical transition can be described in a static framework by introducing 
a replicated version of the system [9, 16]. In fact, the replica method allows exactly to compute averages of the form 
(•)^ (•)in' that enter in Eq. (31), in a static framework without the need of solving the dynamics. For every particle 
we introduce m — 1 additional particles identical to the first one. In this way we obtain m copies of the original 
system, labeled by a = 1, . . . , m. The interaction potential between two particles belonging to replicas a, b is Uafc(f ). 
We set Vaaif) = vir), the original potential, and we fix Vabir) for a ^ 6 to be an attractive potential that constrains 
the replicas to be in the same metastable state. 

Let us now define our basic fields that describe the one and two point density functions 

AT 

Paix) =y^J{x-xt) , 
Pab i^^y) Pa{x)pb{y) - Pa{x)SabS{x - y) . 

To detect the dynamical transition one has to study the two point correlation functions when Vab(r) — >■ for a ^ b, 
and in the limit m — >• 1 which reproduces the original model [9, 16]. We denote by the equilibrium average for 
the replicated system under the conditions stated above. The crucial observation is that in the limit Vab{T) ~> 0, 
all replicas fall in the same state but are otherwise uncorrelated inside the state. This leads to the following rule to 
compute the average one should 



' " ' (33) 



• replace (•)^ = (•)^, 

• factorize the averages i*)^^ when they involve different replicas, and 

• remove the replica indexes. 

For instance, for any spatial argument, and for a 7^ 6, we have that following the prescription above 



{PaPb)r = {PaMn, = (Pa),^ (pb) n, = (P) m (P) m (34) 

which is exactly the kind of average we want to compute. Similarly, assuming that different letters denote different 
values of the indexes: 



{PaPaPb)r = (PaPaPb),,, = {PaPa)^ (P&>m = iPP)m (P) m ' ,orN 

(35) 

{PaPbPc)r = {PaPbPc)^ = {Mm (Mm (Mm = (p) m (p) m (p) m 

Let us introduce a space-dependent order parameter 

qab{r) = I dxfix)p';^,\r-x/2,r + x/2) , (36) 



and the two-replica correlation function 

Cab{r)^{qab{r)),-p^^ f dx f{x)[{pl^^ {r - x/2,r + x/2)), - {pa{r ~ x/2)),{pb{r + x/2)),] , (37) 
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where f{x) is once again an arbitrary short ranged function. We are interested in these functions for a ^ b. Using 
the prescriptions above we obtain 



Cab{r) = J dxf{x)[{pa{r~x/2))^{pb{r + x/2))^]-p' . (38) 

At this point the rephca indexes can be dropped because the one-rephca average in a metastable state is the same for 
all replicas, and we get 



Cab{r) = J dx fix) [{p{r - x/2))Mr + xl2))^ ^ p' = E[(C(r, t ^ oo))] - p' , (39) 

which provides the crucial identification between replicas and the long time limit of dynamics in a metastable state. 

Similar mappings can be obtained for four-point correlations. We define the correlation matrix of the order param- 
eter as (for a b and c d): 

GiiL(^) = (9a6(r)gcd(0)>, - (qabir)), (qcM), , (40) 

where the superscript / is useful to keep in mind that we performed a smoothing through the function f{x). Performing 
similar manipulations as for the two-point functions, we have 



GiiUW = J ^'''^y fi^)f(y) [(/5a - x/2)pb(r + x/2)p,i-y/2)p,{ymU - E[(C(r, t ^ ^))]E[(C(0, t ^ oo))] (41) 

and in the first term the average (•),-,-j can be factorizcd over different indexes. Comparing this with Eq. (32) we 
obtain: 

(f) (f) (^^2) 

G,.(r,i->cx3)=Gi{|,,(r)-Gi{|,,(r) . 

D. The replicated free energy 

We now discuss how replica correlation functions can be computed. We introduce some standard notations of liquid 
theory [36] and we adapt them to the replicated system. 

Let us start with the grand canonical partition function for a _D-dimcnsional fluid with pairwise additive potential 
V, chemical potential p, and under an external field ^E*. The logarithm of the partition function reads: 

W[iy,w] =\nZ[i^,w] =\n^ ^1 J^dx^ exp f dxdy p^^\x,y)w{x,y) + j dx i^{x)p{x)j (43) 

where we have used the following definitions for the fields 

N 

(44) 

P^'^Hx^y) ^^^6{x - x^)6{y - Xj) = p{x)p{y) - p{x)5{x - y) 

i=l j^i 

and the microscopic details of the system are encoded in 

Ka;)=/3/i-/3*(x) , ^^^^ 
w{x,y) = ~f3v{x,y) . 

To study the glassy phase we will follow the method introduced in [9, 16]. We replicate the system introducing 
other m — 1 copies of this original fluid, with interaction between copy a and copy b denoted by Wab = —Pvab, so that 
the logarithm of the replicated partition function is given by 

W[{iya}, {wab}] = In Z[{iya}, {Wab}] = (46) 

oo -. „ / rn N \ /-. l,m „ m „ 

= E / n n^^^^"^ 2 E / d^dy p'^2{x.y)wab{x,y) + dx ua{x)pa{x) I , (47) 
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where the definition of the fields must be modified in order to take into account difl^erent rcpficas 

N 

Pa{x) = ^S{x - XiJ , 



4=1 



(48) 



(2) 

In the fohowing, to hghten the notations, we will sometimes (when this leads to no ambiguity): 

1. use shorthand notations for the spatial positions, e.g. /(I) f{xi), /(I, 2) f{xi,X2)', 

2. similarly, use Jdxi — > J-^; 

3. drop the replica and space indexes and simply denote p — )• {/5o(l)}, P^^' {plibi^^ 2)}; 

and similarly for similar or more complex quantities. 

To study the glassy phase we need to know the correlation functions of the two fields p and /S*-^-* . Let us underline 
that in this scheme, the details of the (attractive) interaction between different replicas is not important because in 
the end we will send it to zero. In fact if there is a glassy phase below at a certain dynamical temperature Td, this 
infinitesimal attractive potential is enough to let all the m replicas fall down in the same state. Hence, it is very 
convenient to perform a double Legendre transform and write the free energy as a function of the averages of p and 
p^'^\ which we denote p and p'-^-' [19, 37, 38]. In this way, we can take directly the limit where there is no interaction 
between different replicas and look for a solution where the replicas remain correlated in this limit. We obtain [9]: 



r[p,p(')] = E / pa(i)^:(i) + ^E / pil\h2)w:,{i,2)-wy,w*] 

a •' a,b -^l'^ 



where v* and w* are the solution of the two equations 



pa{l) and 



5W[y, w] 



^Pil^(l,2). 



Morita and Hiroikc [37] showed that this double Legendre transform can be written as 

r[p, p(2)] = rid[p, p(2)] + r,i„g[p, p^^'] + r2Pi[p, p^^^ , 



(49) 



(50) 



(51) 



where 

rid[p,p(^' 



ring 



[P,P 



(2)1 



J2 I Pa{l)[\r.pa{l)^l] + W I 

1 (-1)" 



Pi'^(l,2)ln 



Pa,b 



1,2) 



pi'^(l,2) + p„(l)pfc(2) 



Pa(l)pb(2) 



E 



n>3 

and we have introduced 



(52) 



{1,2) ■■■ pa„{n)ha„atin,l) 



habil,2) 



pS(1:2) 

Pa{l)pb{2) 



(53) 



The term r2pi is the sum of all two-particle irreducible diagrams that are defined precisely in [19, 37, 38]. We do not 
give more details because this term will be mostly neglected when we will perform concrete numerical computations. 
However, the formalism we develop below holds in full generality so one could include 2PI diagrams in future works^. 

It will be useful in the following to define the direct correlation function Cab(l,2) through a replicated version of 
the Ornstein-Zernike equation: 



hab{l,2)=Cab{l,2)+y] [ /la,(l,3)pc(3)c,b(3,2), 
„ "'3 



(54) 



^ For instance, this can be done through a systematic expansion in powers of the off-diagonal term of p^^^ . The price to pay is that the 
result depends on many-body correlations of the non-replicated liquid [39] . 
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whose solution can be written through a series expansion in the foUowing way 



:afa(l,2) =;^(-l)" J2 I /iaai(l,3K,(3)V„ 
,1^1 a2,...,o„_i "'3---"-l 



, (3, 4) • • • pa„_, {n - l)/ia„_ifa(n -1,2). (55) 



The free-energy is computed by evaluating the functional F at the physical correlator p^^(l,2) which solves the 
following equation [19, 37, 38]: 



(56) 



Pab 



The density field can be determined by a similar equations as a function of the chemical potential, however here we are 
interested in a solution with Pa(l) = P, hence we can directly fix the density in this way. Moreover, we are interested 
in a solution for the two point function which has eventually (below the dynamical temperature Td) a IRSB structure 



P(^)a6(l,2) 



W5(1,2) + (1-Mp'5(1,2) 



(57) 



In particular in the high temperature phase we expect that the off-diagonal part of this solution, namely g, is trivial 
(it corresponds to uncorrclated replicas, hence it is identically equal to 1) while below the dynamical temperature 
Td we have a non trivial solution. Note that the glass transition can be crossed either by lowering temperature or 
by increasing density, the second strategy being the only possible one for hard spheres like systems. In the following 
general discussion we will typically refer to lowering temperature: but all of our results apply to any other path in 
the phase diagram that crosses the glass transition line. We will indeed present concrete numerical calculations both 
in temperature and density. 



IV. LANDAU EXPANSION OF THE FREE ENERGY AROUND THE GLASSY SOLUTION 



In this section, we will show how one can construct a Landau expansion of the free energy. There are two main 
differences with respect to the simple ferromagnetic example of Sec. II. First of all, even for a homogeneous system, 
the order parameter is g{xi — x-i) and it keeps a non trivial dependence on space. Second, at the mean field level 
the glass transition is a random first order transition: even if it is a second order transition from a thermodynamic 
point of view, the order parameter has a finite jump at the critical (dynamical) temperature Td- This implies that we 
cannot approach smoothly the glass phase from the liquid one. In the following, we will assume that we are in the 
glass phase at T < Td, where Tj is non-trivial, and study how the limit of e = — T — > is approached from positive 
e. 

In Sec. IV A we introduce an appropriate scalar order parameter and perform a Landau expansion of the free energy 
for small deviations of the order parameter around the critical point. In Sec. IV B we show how this expansion can be 
used to compute the MCT exponents, in particular the exponent A, following [14]. In Sec. IV C we perform a more 
detailed study of the mass matrix, i.e. the quadratic term of the expansion. In Sec. IV D we use this to show that 
the value of the MCT critical exponents do not depend on the details of the definition of the scalar order parameter; 
as a side product we obtain a much simpler expression for these exponents. 



A. Free energy for a uniform field 



We want to define the free energy as a function of the order parameter Ca6(r), defined in Eq. (37), in the case in 
which it is uniform, i.e. independent of r. The way we can produce this quantity is just by maximizing the free energy 
with respect to p^^^ under the constraint that Cab is given by its definition, which is enforced through a Lagrange 
multiplier Eat- 



^[Cab] = max 



r[p,P^'^]-^eab ( VCab- 
a^b 



f{l,2)[pll\l,2)-pail)pbi2y_ 



1.2 



(58) 



We now that in absence of the constraint, Sab ~ 0, the free energy F[p, p^^^] is maximum for p^-^' = Eq. (57), 
which corresponds to some value Cab of the order parameter. We want to expand around this reference solution. We 



caU Ap 



(2) 



P 



(2) 



and ACab — Cab — Cab thc deviations from this reference solution. 
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We use from now on a lighter notation in which a global index is used instead of replica and space indices, 
A = {a, &, 1,2}, B = {c, d, 3,4}; moreover we use the Einstein's convention where repeated indices are implicitly 
summed. Expanding r[/ci, p'^'] around the reference solution we obtain (recall that p — Pa(l) = Ph(2) is a fixed 
constant) 



AT[ACab] = max 



1 (,) (2) 1 S^r[p,pi^)] . J2) (2) (2) 

^Pa ^Pb + 6 c (2). (2). (2)^^a ^Pb ^Pc 



^ea6(^yACa6-^^/(l,2)Apf,'(l,2)) I , 



(59) 



where all the derivatives must be computed in the reference solution. Defining 



Mab 
Labc — 



X (2)x (2) ' 

SPa ^Pb 

e (2)r (2)r (2) ' 

^Pa^Pb^Pc 



(60) 



and ea = £abf{^, 2) (where we define Saa = 0), the derivative with respect to Ap*^^^ leads to the following equation: 



= MAB^pf + -LABcAp'^'Ap'-c + £A 



1 



.(2) A J2) 



which can be inverted perturbatively and gives 



Plugging this in the free energy we obtain 



Ar[ACab] = max 



^^^A-^'^AB^S — -^LABcMj^\,Mg^,M(j^,£A'£B'£C' — ^ ^ £ab ACafc 



a^b 



(61) 



(62) 



(63) 



At this point it is convenient to recall that £a = £abf{^,2). Using the shorthand notation {M ^L)abc' 
^AA'^^BB'^^cc'^A'B'C and introducing 



I,- - ,6 



VCab,cd,ef 

we can rewrite Eq. (63) as 

Ar[ACa6] Fmax 



/(l:2)M-yi,2,3,4)/(3,4) , 

/(I, 2)/(3, 4)/(5, &){M-Hr^M~^L)abM^,f{l, 2; 3, 4; 5, 6) 



(64) 



"^^abMab.cd^cd ~ 7:^ab.cd,ef£ab£cdSef — ^ EabACafi 



We now take the derivative with respect to e^b a-nd we obtain 



which is inverted as 



ACa& = —-^ab.cd^cd — ■^^ab,cd,ef£cd£ef , 



Eab — —Mab,cdA-Ccd — -^M.ab,cdC-cd,ef,ghM.ef,e'f'A.Ce'f'M.gh,g'h'A.Cg'h' , 



(65) 



(66) 



(67) 



and plugging this in Eq. (65) we finally obtain the desired third order Landau expansion of the free energy, which is 
the analog of Eq. (5) for the glass transition: 



Ar[ACa6] = V S^^ACabMab^cdACcd + ^Wab,cd.ef ACabACcdACef 
yVab,cd,ef = Mab,a'b'Mcd,c'd'Mef,e'f'jOa'b',c'd',e'f' ■ 



(68) 
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B. Computation of A 



At this point we are equipped with aU the ingredients to give an exphcit expression for the parameter exponent A 
of MCT. It has been shown in [14] that in mean field disordered systems this parameter can be computed in a purely 
static framework. The argument is the following: the dynamics of such systems can be studied in great detail using 
the Martin-Siggia-Rosc [40] formalism; in particular, by going to the supersymmetric representation for the action of 
the Langevin process describing such dynamics [41], one can obtain a dynamical Gibbs free energy which has exactly 
the same form (apart from a kinetic term containing the derivatives with respect to time which can be neglected if 
we want to study the long-time behavior) of the one in (68) where replica indices are replaced by supertimes. From 
the dynamical action it is straightforward to see that the exponent parameter A is given by the ratio between two of 
the cubic terms of the dynamical Gibbs free energy. But now the key point is that this two coefficients are related to 
six-points dynamical correlation functions whose long time behavior can be studied using replicas following the same 
analysis of Sec. Ill C. In particular it can be shown that for times such that the two point correlation function is very 
close to its plateau value, the values of these dynamical six-points correlation functions can be computed in a static 
framework just by using the replicated Gibbs free energy of the form (68). This implies directly that the exponent 
parameter can be computed from the statics. This argument has been discussed in full detail in [42]. 

The cubic coefficients of the replicated Gibbs free energy which are relevant for the computation of the A are the 
following ones 

-^Tr(AC)3-!^5:AC3,. (69) 



and the exponent parameter is given by 



^ • (70) 

Wi 



In computing the derivatives in (60) one has to use the explicit form for the order parameter pf^) that maximizes 
the free energy, given in Eq. (57), which is replica symmetric (corresponding in our formalism to a IRSB solution). 
Exploiting this symmetry, in [43] it has been shown that for a replica symmetric saddle point the two coefficients wi 
and 1X12 can be written in the following form 

wi^Wi- 3>V,5 + 3W7 - 8W8 , 

1 3 (71) 

'W2 = - 3W3 + -Wi + 3W5 + 2W6 - 6W7 + 2W8 , 

where 

Wl = Wab.bcca , 
W2 = Wab,ab,ab , 
W3 = Wab,ab,ac , 

W4 = yVabMb,cd , , , 

W5 = Wab^acM , ^ ' 

We = Wab^acad , 
W7 Wac,bc,de , 
Ws = Wab,cd,ef ■ 



Therefore the final expression for A is 

?«2 _ ^W2 - 3W3 + |>V4 + 3>V5 + 2W6 - 6W7 + 2^8 
Wi - 3W5 + 3W7 - 8W8 ■ 

Note that this expression contains implicitly a dependence on the function /(x) that we have chosen to define the 
order parameter in Eq. (37). Moreover, note that the usefulness of this relation is very hard to prove at this stage 
because to compute the cubic coefficients W we have to perform a complex operator inversion to obtain the M., see 
Eq. (64), followed by a convolution with the £, see Eq. (68). This would be an extremely hard numerical calculation. 
In fact, all these problems can be solved noting that the mass matrix of the free energy develops a zero mode at the 
dynamical transition, and the presence of this zero mode simplifies a lot the computation of A, as we show next. 
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C. The mass matrix and the zero mode 



Recall that here we are looking to a homogeneous system. Hence, in Eq. (57) we have g{xi,X2) = g{xi — X2). 
Keeping this in mind, in the following we will use equivalently the notations g{xi^X2), 5(1,2), and g{x). As we 
already explained, at the mean field level the solution of the saddle point Eq. (57) is g{x) = 1 in the liquid phase 
while ^(x) is a non-trivial function in the glass phase. The appearance of an off-diagonal solution below the dynamical 
transition point is discontinuous and can be regarded as a bifurcation-like phenomenon at some critical value of the 
control parameters (temperature or density). Because the transition is discontinuous, the bifurcation does not happen 
around the liquid solution g{x) = 1, but around a non trivial value of g{x) corresponding to the solution at the critical 
point. Therefore, if we approach the transition point from the liquid phase and look at the behavior of the free energy 
around the liquid solution, nothing special will happen. For this reason we are forced to take the control parameter 
such that we are in the glass phase (for example we put T < Td) and approach the transition from this phase. Let us 
call e = Td~ T the distance from the critical point. Then for e ^ 0+ 



5(1, 2; e) ^ 5(1, 2; 0) + 2 At fco(l, 2) + 0(e) 



(74) 



Here, the function fcg is normalized by ^ /i 2 ^o{^, 2)^ ~ J da;fco(a;)^ = 1, which defines implicitly the constant k. 
Therefore 



^a6(l,2) 



^il-Sab)Kko{l,2) + Oil) 



de y/e 

We consider the saddle point condition (56) with a =^ b and Wa^t = 0. Then 



(75) 



E 



= E 

Pab ^=^<^ 



3,4 <Spf,^(l,2)5pr;(3,4) 



(2) 



rfp,d(3,4) 



de 



Pab 



hA <5pf,)(l,2)5pi?)(3,4) 
Recall the definition of the "mass operator" in Eq. (60) 



dp,,(3,4) , S'T[p,p' 



(2) 



de 



5p^\h2)Se 



Pab 



(76) 



Pab 



b:cd 



(1,2; 3,4) 



<52r[p,p(2)] 



^Pi'^(l:2)<5pf,^(3,4) 



(77) 



Pab 



for a ^ b and d. Then we can write Eq. (76) as 



3.4 



AW,(1,2;3,4)^^^ + /C(1,2) 
de 



(78) 



where /C(l, 2) is finite at the critical point and does not depend on a 7^ 6 because of the symmetry of the saddle point. 
Recalling from Eq. (75) that the derivative of the saddle point solution is divergent when we approach the dynamical 
point from below, we conclude that the mass operator should develop a zero mode at the transiton, in such a way 
that the divergent part of Eq. (78) is cancelled. In other words, we should have that 



E 



Afab;e<i(l, 2; 3, 4)fco(3, 4) = M\/^fco(l, 2) 



(79) 



3,4 



Because of the replica symmetry of the saddle point solution, the most general form for this mass matrix is [43] 



M,b;crf(l,2;3,4)-Mi(l,2;3,4) 



SacSbd + SadSbc 



1/2(1, 2; 3, 4) 



Sac + Sad + She + hd 



il/3(l,2;3,4) (80) 



Because such matrices make a closed algebra, the inverse of the mass operator must have the same replica structure: 



Ga6;cd(l,2;3,4) = Gi(l,2;3,4) 



SacSbd + Sad^bc 



G2(l,2;3,4) 



Sac + Sad + Sbc + Sbd 



G3(l,2;3,4) 
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The equation for G is 



! AW(1,2;3,4)G.«(M;5,6) ^ ^.-^^Z + ^(1.5)^(2. «) + ^(1,6)^(2.5) ^^^^ 

and, in the repUca Umit m — > 1, the solution of this equation is given by^ 

Gi = 

G2 = -2[2Mi - M^r^ ® Af2 ® (83) 
Gg = Aff ^ ® {M2 ® [2A//i - A'hr^ ® Af2 - M3} ® Mf ^ . 

We know from Eq. (75) that the zero mode is independent of the off-diagonal replica indexes. This implies for a 
matrix of the type (80) that in the replica limit m — > 1 



lim V / AWd(l,2;3,4)fco(3,4) 



lim / [m(TO-l)A/3(l,2;3,4) + Mi(l,2;3,4) + (m-l)Af2(l,2;3,4)]fco(3,4) (84) 



3,4 



A/i(l,2;3,4)fco(3,4) 

This implies that for m ~ 1 among all the components of the mass operator, fco is a zero mode (with eigenvalue 
proportional to y^) only for the operator defined by the kernel Afi(l, 2; 3, 4), while Af2 and M-s do not need to have 
any zero mode at the transition. In formulae, from Eq. (79), 



A/i(l,2;3,4)A;o(3,4) =/xViA:o(l,2) , (85) 
or equivalently (recall that ^ fco(l, 2)^ = V): 

Mri(l,2;3,4) = -^^fco(l,2)A;o(3,4) + 0(l) . (86) 

This observation is very important because it shows that the operators defined by the kernels Gi and G2 have a single 
pole (a divergent eigenvalue) while the operator G3 has a double pole. This is a very straightforward generalization 
of results obtained in [13] to the case in which the system has a non trivial spatial structure. Let us now rewrite 
Eq. (83) as 

G2 = 02®Arr' with O2 = -2[2A/i - M2]^i ® M2 , . 

Gz = M^^ ®Oz®M:[^ with O3 [A/2 ® [2Afi - A/2]"^ ® Af2 - Af3] ^ ^ 

and use this to conclude that the most divergent contributions to these kernel operators are given by 

Gi(l, 2; 3, 4) ~ ^fco(l, 2)fco(3, 4) + 0(1) 

G2(l, 2; 3, 4) ~ ^fc2(l, 2)fco(3, 4) + 0(1) with ^2(1, 2) = ^ 02(1, 2; 3, 4)fco(3, 4) 
G3(l, 2; 3, 4) ~ p^^fco(l, 2)fco(3, 4) + o(-i=) with ^ ^ ^ ^ ^ fco(l, 2)03(1, 2; 3, 4)fco(3, 4) 



Note that in G3 a term proportional to 1/ ye appears and it depends on the exicited states of the kernel operator Afi . 
However we will see that the contribution of the zero mode is enough for the computation of A. 



Wc denote with ® the convohition in the space variables 

a(gjb= f a{l,2,; 3,4)6(3,4; 5, 6) 
J'iA 
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D. Analysis of the cubic terms and computation of A 

At this point we are equipped to extract the divergent part of the cumulants wi and W2 and obtain a simple and 
universal expression of A. Let us start with the generic expression (68) for the cubic coefficients of the free energy as 
a function of the a uniform order parameter. Using Eq. (64) we have 



Wab,cd,ef ^ Mab:a'b'Mcd:c'd'Mef;e'r^ f /(I, 2)M-,\.^„,„ (1, 2; l', 2') X 

Jl,...,6;l',...,6' 



X /(3, 4)A/-],.^„^„ (3, 4; 3', 4')/(5, 6)M-},.^„f„ (5, 6; 5', 6')La"b";c"d";e" /" (l', 2'; 3', 4'; 5', 6') (89) 

^ ^ Aab;a"b" (l', 2')A,d-c"d" (3', 4') Ae/;e",/" (5', 6')2va"6";c"d";e" /"(!', 2'; 3', 4'; 5', 6') , 



where we introduced the quantity 

Aab;a"b" (l', 2') = Mab;a'b' f /(I, 2)M-,l^„,„ (1, 2; l', 2') 



(90) 

= Ai(l', 2')^-^-"^'-'^" + A2(l', 2^)W+W+'?a.b"+'^M" ^ ^^(^,^ 2') . 

Clearly the replica structure of the matrix Aiab-cd is inherited from the structure of the mass matrix Mab-cdi^, 2; 3, 4) 
and its inverse. From Eq. (64) we have 



^ f SacSbd + SadSbc\ , ^ f Sac + Sad + She + Sbd\ , ^ 
^ab.cd = ^1 ( 2 ) + ^2 I \+Qz , 

.. ,. I SacAd + SadSbc \ , , , f Sac + Sad + Sbc + Sbd\ , 
Mab.cd = Ml \ +M2\ + M3 



2 

where we have defined the following /-dependent quantities 



(91) 



V Ji 

and where using the same algebraic structure of the replica sector 

Mi = ^ , 



G. = ^ I /(l,2)G,(l,2;3,4)/(3,4) , (92) 

2,3,4 



2^2 

^ ' gr(2g.-g2) ' (9^) 



M2 



^^(2^1- ^2) Gi ■ 

Next we have to analyze the matrix A that appears in Eq. (89); performing the matrix multiplications we obtain 

Ai(l',2') = /*^Gi , 
yi 

01(201-02) ' (94) 

0!{20i~02) ' 

where we used the shorthand notation (/*Gi)(l', 2') — Ji 2 2)Gi(l, 2; 1', 2'). Now we can use the critical behaviour 
of Gi in Eq. (88). We see that divergent terms of the same order will appear at the numerator and denominator 
of the above expressions, so we can extract the finite part of A^ in the limit e — > 0"''. Thanks to some non-trivial 
cancellations the result is, defining f * kg ~ V^^ ^ /(li 2)fco(l, 2) — J dxf{x)kQ{x): 

J * ko 

A2(l,2) = 0, (9^) 
A3(l,2) = , 
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and plugging this in Eq. (89) the final expression for the cubic coefi'cicnts is given by 

mb-.cd-.ef - \ ,3 -^ / fco(l, 2)fco(3, 4)fco(5, 6)ia&;cd;e/(l, 2; 3, 4; 5, 6) . (96) 

(J * Koj V Jl,2:3,4;5,6 

Using Eq. (71) we obtain the two relevant coefficients for the computation of A as: 

^2 = 77 rr^TP f ko(l,2)ko{3,4:)ko{5,6)(-Lab,ab.ab~'iLab.ab,ac+-zLab,ab,cd + 3Lab,ac,bd + '^^ab.ac,ad 

(/★fco)'^ 1/ 7l 2;3,4;5,6 V2 ■ 2 

~ (^Lac,bc,de + '^Lab.cd,ef (1,2; 3, 4; 5, 6) (97) 

1 If 

wi = -7-7 — rTaTJ / ^o(l, 2)fco(3,4)A:o(5, 6){Lab.bc,ca - 3Lab,ac,bd + iLac,bc.de - Lab,cd.ef){i, 2; 3,4; 5, 6) . 

U * '^O) V 7l,2;3,4;5,6 

Note that because A is the ratio between W2 and wi, the dependence on / disappears and we get a consistent /- 
independent expression for the exponent parameter. Moreover, the expression above are much simpler to compute 
because they simply require a convolution of the functions L with the zero mode fcg, without any need for operator 
inversion. 



V. GRADIENT EXPANSION 



We want now to study the critical behavior of Gth{r) and 64 (r). These quantities are related by Eq. (42) to the 
inverse of the mass operator G'af,;cd(l, 2; 3, 4), Eqs. (80) and (81). In principle we should invert the mass operator 

( f) 

Mab;cd{i, 2; 3, 4) and perform a convolution with the smoothing function f{x) to obtain the replica correlation G)j^.^^(r) 

that enters in Eq. (42). Although certainly well defined, this procedure would be numerically quite heavy even in the 

simplest approximation for the mass matrix (e.g. the HNC approximation we discuss below). 

In this section we want to show that, if we are only interested in the long distance behaviour of G)j(j^^(r), we 

can take advantage from the existence of the zero mode discussed in Sec. IV C to obtain a universal critical form of 
f f) 

G)j^.^^(r). In particular we will show that the function f{x) only enters in the prefactor, while the correlation length 
is independent of f{x). For this we need to set up a gradient expansion where we consider a field Ap^jj,^ (1, 2) which 

f2l f2) 

is almost uniform. Recall that in the uniform case Ap^^^ (xi, 0:2) = Ap^^ (xi — X2)- Here we consider a non-uniform 
field, that depends also on the variable (xi + X2)/2, but we consider that the dependence on this variable is weak. 



A. Fourier transform 



It is convenient to separate the dependence on the translationally invariant variable xi~X2 from the slow dependence 
on space (xi + x-i)!^ by introducing a Fourier transform 

Apl''(p,g) = / dxid.T2e^^'(^)+^«(— ^)Apf,)(xi,X2) , 

Apil^(-i,X2) = /^e-(-^-)--K-)A,i^)(,,,) . 

Here p is the momentum coupled to the slow spatial variation, while q is the momentum coupled to local displacement. 
Hence we will be interested in the limit of small p. Plugging this in the quadratic part of the expansion of the free 
energy around the saddle point solution, we obtain 



A2r=^ E / Api'^(l,2)AW(l,2,3,4)Apf;(3,4) 

^ a^b c#d-^l'2,3,4 



1 /^i^^A,i^)(p,,)M..,.(-p,-^;-,,-fc)Apf,)(p,fc) 



(99) 



where 



Mab-cd{p,P:q.k) = / e'H^)+^9(-i-=)+'p(^)+''=(-3--4)^^^^.^^(1^2;3,4) (100) 

Jl,2,3,4 
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Because of the translational invariance of the saddle point solution, the mass matrix is also translational invariant, so 
that we can make a change of variable to X = {xi + X2 + + X4)/A and Ui = Xi — with X)^=i ~ 0' ^^'^ then 

M(l, 2; 3, 4) does not depend on X. Calling Vu — duidM2du3du4i5 Si=i '^^ S^t 



(101) 



We want to study the correlation function 



Vu e 



(102) 



when the two points ui, it2 are far away from M3, U4. This implies that as already announced we need to study 
the limit p — > and we have to develop the mass matrix around this limit. The correlation function above, at the 
Gaussian level, is given by the inverse of the mass matrix, i.e. by transforming Eq. (82) to Fourier space we get 



^ n \^'^) n 



(103) 



Therefore it has the form, akin to Eq. (83), 



(P) 



M- 



(p) 



Gi^^ = -2 



(g) M^^^ <Si 



Ml 



(p) 



(104) 



a 



(p) 



M- 



(p) 



(g) i M. 



(p) 



,f{p) 



The convolution and inversion operations appearing in Eq. (104) are defined by 



(105) 



Let us now consider the four point function Gj^.^^{r — r') defined in Eq. (40), which we can write as 



Gl]cd{r ~r')= dxdy f{x)f{y) Gabbed (r + |, 



X , ^ y , y 

r ;r H — ,r 



(106) 



We have 



{27T)^S{p + p)G'-/l,{p)^ I drdr'e^P^+'P^' I dxdy f{x)f{y) G^b^ (r + |, r - |; r' + |, r' - | 

Sj/(-g)./(-fc) / drdr'dxdye'P-+'P'''+''>^+'''y Gab,cd + 



dg dk 
dq dk 



X X , y , y 
2'" +2'^^ -2 



f(-q)f{-k) [ e'P(^^)+p(^^)+'"^(^i-^^)+''=(^^-^*)G'af,;cd(l,2;3,4) 



1,2,3,4 



{2n)^S{p + p) I ^^f[-q)f{-k)G%MM ■ 



(107) 
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Hence we finally obtain 



Using Eq. (42) we find that 
dq dk 



(108) 



Gthip) = 



f{-q)f{-k) G^:l,iq,k)-GZ'{l^k) 



(p) 



(27r) 



2D 



lG[^\q,k) + lG^^\q,k) 



(109) 

We can now perform some formal manipulations on the operator, dropping the indexes for convenience. Using 
Eq. (104) 



2 4 All 
_ 1 
" Ml 
1 

~ Ml 
1 



1 



M2 


- 2Mi 




1 


Ah 


- 2Mi 




1 


M2 


- 2Mi 




1 



M2- 1 

M2- 

2Mi 



1 



2Mi 
1 



M2 - 2Mi 
1 



{AI2 - 2 Ah) 



2Mi 



(110) 



2A/i 



Ml M2 - 2A/i 

so that the critical behavior is driven by the dominant divergent contribution to the kernel operator Mi: 

-1 



Ai^ 



(p) 



(111) 



B. Spectrum of the mass matrix 

Let us now study the kernel operator Mi in more detail. We have seen in the previous sections that it develops 
a zero mode when approaching the transition from below. Moreover we know that the zero mode is translationally 
invariant. Let us now look at the spectrum of this operator. The crucial observation is that because of translational 
invariance, the mass matrix is proportional to S{p + p), so we can diagonalize the kernel AI^^\q,k) at fixed external 
momentum p. We can thus write a spectral decomposition 

A4^\q, k) ^ Ao(p2)4'^(g)4'^(fc) + Up')i^i''Hq)i'^lHk) (112) 

a>l 

where we put the lowest mode in evidence for future convenience and 

.M(^)(q,fc)^(f)(-fc) = A„(p2)^(p)(g) , 



(27r)^ , , 

j^i^i'Hk)^l^}i^k)^6^,^, . 

In what follows we will suppose that there is a persistent mass gap between the zeroth eigenvalue and the first excited 
ones, even at the transition. This means that even if the zeroth order eigenvalue goes to zero at the transition, we 
are assuming that all the excited ones remain finite. Let us now consider again Eq. (85) for the zero mode 

Mi{xi , X2, X3, X4)ko{x3 - 2:4) = (Ai\/e + 0{e))ko{xi - X2) ; (114) 

3,4 



in Fourier space it is given by 

dfc 



■A4''-°\q,-k)ko{k) = {^l^+0{e))koiq) . (115) 



(27r)^- 
Therefore 

i'i'^^^q) = fco(g) , Xoip ^0)^fiVe + 0{e) , (116) 
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Let us now define (for i = 1, 2, 3) 



dq dk 
(27r)2-D 



ko{-q)M^P-''\q,k)koi-k) 



in such a way that 



(117) 



(118) 



Hence Ao(p^) is small close to the transition and for small momentum p ~ 0, so that the leading term in Aff ^ is the 
zero mode: 



M. 



(p) 



{q,k) 



Q> 1 



1 



■4^Hg)4^Hft) (119) 



We now need to compute the coefficients /i and a. The first one is given by the definition 



dmi 

^ = hm = lim 



d /"d^qd^fc 



d^e d^ 7 (27r) 



fco(-<7)A/^^-"'(<7,fc)fco(-fc) , 



(120) 



and the other one can be evaluated using perturbation theory. In fact, let us develop the operator Mi around p — 
We have that 



M'f\q,k) = Ml''=">{q,k) ^M^'>iq,k) 



d 



Ap)/ 



(121) 



p=0 



and we know that ko{q) is the fundamental state of the unperturbed operator M[^~^\q,k). Now let us treat the 
second term as a small perturbation because we are near p = 0. The shift of the eigenvalue of the ground state is 
given by 



Ao(p^) - mi 



koi-k) , 



p=0 



so that 



(122) 



(123) 



p=0 



We conclude from this analysis, and in particular from Eq. (118) and (119), that close to the transition and for small 



M 



ip) 



1 



;ko{q)ko{k) 



Inserting this in Eq. (Ill) we get, with the same definition of the ★ product as before: 



In this way it is clear that the correlation length is given by 



dq 



(27r) 



D 



.f{-q)ko{q) 



and 



Go — 



-1/4 



/i 



(124) 



(125) 



(126) 



(127) 



This completes the analysis of the quadratic part of the action. 
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VI. GINZBURG CRITERION 



111 the previous section we have investigated the quadratic and cubic terms in an expansion of the free energy 
around the critical value of the order parameter. This corresponds to the Landau expansion described in Sec. II A 
and Sec. II B. We now follow the discussion of Sec. II C: we assume that the free energy F has been obtained by 
truncating in some way the high temperature expansion (we will come back on this point in Sec. VII), and we use 
this mean field free energy as a bare action to perform a loop computation and obtain a Ginzburg criterion. 

The quadratic part of the free energy provides the Gaussian part of the bare action. What can be expected from 
the discussion of Sec. II D is that the Gaussian approximation is valid (at least qualitatively) everywhere in the (T, p) 
plane if the dimension is greater than the upper critical dimensions; otherwise the Gaussian approximation is valid 
only up to a certain distance from the critical line in the (T, p) plane. As discussed in Sec. II D, this statement can be 
made more precise by performing a one loop expansion for the correlation function of the order parameter and then 
look at which temperature or density the first loop correction is of the same order of the Gaussian approximation for 
that correlation. The Gaussian bare propagator is G''^^.^^{q,k) given in Eq. (104). To compute the loop corrections 
to this propagator, one should start a very difficult computation because the diagrammatic arising from the replica 
field theory is complicated by the presence of the replica indices as well as the "internal" wavevectors (7, k. 



A. Projection on the zero mode 

The calculation can be greatly simplified if we focus only on the critical behavior of the correlation functions. In 
particular, in Sec. VB we have shown that the critical part of the propagator G^ab-cA^^^) entirely dominated by 

the presence of the zero mode of the kernel operator M^\q,k). For example, according to Eqs. (Ill) and (124), 
the singular part of the thermal correlation function Gth{p) is given by the inverse of the zero mode of this operator. 
Clearly, the full thermal correlation function contains also the non critical contribution coming from the excited 
states but here we are not interested in this contributions. Instead, we would like to focus on the critical part of 
the observables. To do this we can suppose that the eigenvalues corresponding to excited states arc set to infinity so 
that we have no fiuctuations along the directions orthogonal to the zero mode. This amounts to consider an effective 
low-energy theory where only the fluctuations along the zero mode arc allowed, which means that we choose 



^p^abip^q) 



(l>ab{p)ko{q) ■ 



(128) 



By doing this we obtain a simplified theory that gives us only the critical part of the correlation we want to compute. 
The resulting effective theory is obtained as follows. The quadratic part is obtained by inserting Eq. (128) in Eq. (99). 
The relevant cubic terms according to [13, 43] can be obtained by inserting Eq. (128) in the cubic part of the expansion 
and using similar considerations as in Sec. IV B. The result is 



dp 



^(Atx/i + (Tp2)|0ah(p)|' + ma ^ 



a^b 





2 






^4>ab{p) 


+ m,3 


^(j)ab{p) 


1 


b 




a^b 





J (9^)21^ Mp)(l>bc{p')(f>ca{-P -P) + -yY1 J (2^)2^ (I)ab{p)(f>ab{p)(l>ab{-P - p) 



(129) 



and the coefficients that enters in the above formula are exactly the ones that were computed in the previous sections: 
see Eqs. (117), (120) and (123) for the quadratic part and Eq. (97) (without the factor /★fcg in the denominator) for 
the cubic part. This result is at the basis of our next computation. Let us note that by projecting on the zero mode, 
the momentum structure of the theory has been simplified drastically. 



B. ip'^ theory in random field 

To get the result at the dynamical transition, one should take the m — > 1 limit. In [13] it has been shown that 
the perturbative expansion of the replica field theory of the form (129) can be mapped to the study of a scalar field 
which satisfies a particular cubic stochastic field equation. This mapping has been done using a transformation for 
the fields which is quite close to the Cardy's treatment of the branched polymer problem [44]. However using the 
field theory techniques developed by Parisi and Sourlas [45], it can be shown straightforwardly that the stochastic 
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equation describes also the leading critical behavior of a theory for a scalar field in a cubic potential and interacting 
with a random Gaussian magnetic field. The action of this spinodal field theory in a random field -which we now use 
as a bare action following the discussion of Sec. II C- is given by 



S{ip)^-J dxip{x){-V^ + m^)ip{x) + ^ J dxip%x)+ / dx{ho{x) + Sh{g, A))ip{x) (130) 
where h^^x) is a gaussian random field with variance given by 



hoix)ho{y) ^ AS{x - y) , (131) 

and the coupling and masses are given in terms of the coupling and masses that appear in the replica field theory 
(129) by 

ml = A^\/e/cr , 

9^{wi- W2)/ct'^' , (132) 
A = -(m2 + m^)/i7 . 

In the action (130) we added a counterterm dh, that will be used to avoid that one loop corrections shift the position 
of the critical point. The counterterm 5h can be seen also as a redefinition of the mean value for the random field 
ho. The bare propagator of the theory is, as usual, given by Go(p) = {p^ + ™o)~"^- To compute loop corrections, it is 
quite useful to write down the generating functional for the connected diagrams 

W[J] ==lnZ[J] , (133) 

where the external current J{x) is given by 

J{x) = hoix) +Sh{g,A) . (134) 
Moreover, let us introduce the diagrammatic notation that will be convenient later 

o = • + • (135) 

where 

J(.t) = o ho{x)= • <5/i(.g,A)=: • (136) 

The expansion for W[J] is given (up to second order in g) by the following diagrams 



o 



W[J] 



o 




(137) 



Let's look first to the corrections to the average of {^{x)) . By taking the average over the random field and by marking 
with a blue dot a doubled propagator, these corrections are given by 



o o 



+ 



(138) 



where we have made the hypothesis (to be checked in a self consistent way) that the counterterm Sh is proportional 
to g so that we hav e retained only the diagrams up to second order in g. We want that the critical point is not shifted 
so we want {^{x)) = 0. Therefore at this order in perturbation theory we must have 



o 



(139) 
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from which we see that 5h cx g. Now we look to the correction to the propagator. We have 



o o 



{(p{x)ip{y)) = Go{x ~ y) 



Go{x-v)+- 



where we have used the relation (139). Because we are interested in the most (infrared) divergent diagrams we can 
neglect the second one. However, note that the second diagram is relevant in the ultraviolet regime because it diverges 
sooner as the number of dimensions is increased (actually it is convergent in D < A). Now, let us neglect the second 
diagram and compute the dotted diagram. We have at zero momentum 



2 J i2n)D- 



G{p = 0) = Go(p = 0) + Go(p = 0)^^ / —i-G^{qY . (140) 



From now on we follow closely the procedure of Sec. II D. We invert the above equation to obtain 

Ag2 Aq 1 2 ^9^ 1 



m^ = G- ip = 0)^m,-—J ^^-^—^^m,-—J (141) 

As in the (^^ theory, the mass correction is UV divergent for d > 6 and leads to a non-universal shift of the critical 
point, which corresponds to 



The distance from the critical point is 



^ ' ^ (142) 



(143) 



A52 /-A 1 



and Eq. (141) becomes 



Aff^ r dq f 1 1 



As in Sec. II D we have to impose the condition that the one loop correction is small in such a way that the mean 
field result m|j = t is not altered. A simpler and completely equivalent (apart from a numerical prefactor of order 
one) condition is obtained by considering 

dml 2 J {2n)D [q^ + to|)4 ' ^^^^^ 

and of course we want the second term to be much smaller than one. Written in terms of physical observables (m|j, 
and i), the one loop correction is UV convergent and IR divergent when D < 8 and the renormalized theory exists in 
this case, leading to a universal form of the Ginzburg criterion, which clearly tells us that the upper critical dimension 
is 8. Sending the cutoff to infinity we get^ 

A = i-:^rf4-£)„.r (146, 



rfm|j 4(47r)^/2 V 2 



The computation can be done by using the standard Schwinger representation for the propagator [25] 



1 /"^ 9 2 

= / dae-"(f"+'" ) 

+ Jo 



and by taking derivatives with respect to m 
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and the Ginzburg Criterion amounts to the fohowing condition: 

On the other hand for _D > 8, as before, the correction is UV divergent (and IR convergent) and therefore the precise 
form of the Ginzburg criterion depends on the regularization. Again in the mean field region the correlation length is 
£, — l/mn, see Eq. (126), hence we can write 



VII. EXPLICIT CALCULATIONS IN THE HYPERNETTED CHAIN (HNC) APPROXIMATION 



In order to perform concrete microscopic computations, we now need an explicit expression for the free energy 
T[p, jo'-^-']. As discussed in Sec. II C, there are several possibilities and many of them have been used in liquid theory [36]. 
Most of them can be seen as resummation of the high temperature or low density expansion. A popular one is the 
so-called Hypernetted Chain (HNC) approximation, which amounts to simply neglecting the contribution of 2PI 
diagrams from Eq. (51). This approximation is known to give good qualitative results (although quantitatively not 
very good) for the structure functions around the dynamical transition [9, 46, 47]. Here we report results based on 
this approximation. Better approximations can be obtained by performing a "small cage expansion" [3, 47] or by 
expanding systematically around the HNC equations [39]. Other approximation schemes have been explored [48, 49] 
but they arc not practical for the present purposes. 



A. Derivatives of the free energy in HNC 



Because we need the expansion of the free energy up to the third order we start by computing its derivatives. 
To simplify the notations, here we will perform the calculation without taking into account the symmetries of the 
replica and spatial indices, and we will symmetrize the result at the end. We also use the shorthand notation 
(5a&(l,2) = SabS{xi - X2). 

The first derivative is 



1 



loa 



Pa(l)p6(2) 



+ i[cfc„(2,l)-/ib,(2,l)] 



(149) 



which can be easily proved from Eq. (52) by representing Fi-ing as a sum of diagrams [36]. Using the same diagrammatic 
representation we find 



= [<5,e(l,3) - pc(3)c„,(l,3)] [^,,(4,2) - p,(4)c,,(4,2)] = p,(3)p,(4)ri^j (1, 3)r(^J (4, 2) 



having defined 



(150) 



ri'^(i,2) 



Pa{l) 



(151) 



From this we can compute the second derivative of the HNC free energy, namely its mass matrix: 



(2)1 



'^Plt^(l,2)-5pf;(3,4) 



'5ac(l,3)(5rffc(4,2) 6ad{lA)Stc{2,i) 



Pi'^(1.2) 
1 



Pa{l)pb{2) 



2 Pc(3)pd(4) 



[Sbci2, 3) - pc(3)c,c(2, 3)] [SdaiA, 1) - Pd(4)cda(4, 1)] 



<Sac(l,3)(5dfc(4,2) 6ad{lA)Sbc{2,3) 



Pa{l)Pb{2) 



2ric^(2,3)r^:(4,i). 



(152) 
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Moreover the cubic terms are given by 



'5pi'^(l,2)<5pr,^(3,4)<5p^^/(5,6) 



.(2) 



1 



1 

2p,(3)pe(5)p/(6) 
1 1 fcbc(2,3) 

2pd(4)pe(5)p/(6) <5/ie/(5,6) 
1 
2 

1 r 



[,Sbc(2,3)-pc(3)cbc(2,3)] 



'5cda(4,l) 



6hef{b, 6) 

[,5da(4,l)-pd(4)cda(4,l)] 

1 



-^,,(1, 3),5rffc(4, 2)<5„e(l, 5)<5/,(6, 2)-^ 

bah (1'2)]^ 



rP,)(2,3)r:,:^(4,5)ry,^(6,i) + n:^(2,5)rr/(3,6)r:,:^(4,i 

Using these resuhs, following Sec. IV A, we expand the free energy around the IRSB reference solution, Eq. (57). 
We have Ap'-^^{1,2) = pi^f,^(l,2) -7^^^,(1,2) and AT = T[p,p^^^] - T[p,^]. The derivatives of T and in particular 
the functions F'^^) niust be evaluated on the IRSB solution Eq. (57), hence 



.(2)/ 



.(2) 



(2)/ 



.(2)/ 



.(2). 



(153) 



ri'^(l, 2) = -Sabil, 2) - 6ab[cil,2) - ?(1, 2)] - cil, 2) 



(154) 



and we have introduced for the matrix Caf,(l,2) the same notation as in Eq. (57), where c(l,2) and c(l,2) are 
respectively its diagonal and off-diagonal part. We obtain AF = A^F + A'^F with 



A^F 
A^F 



-y 



1 



- 1 



[Ap,6(l,2)]2 + i Yl I Ap„fc(l,2)rf^)(2,3)Aped(3,4)r 

*a5^b,C5^d -^1.2,3,4 



.9(1,2) 

\ E / Ap„,(l,2)Fi^,)(2,3)Ap,,(3,4)F|,^j(4,5)Ape/(5,6)Ffj(6,l) 



(4,1) 



(155) 



a^b^c^d^e^f 



1,2,3,4,5,6 



B. The HNC mass matrix 



The mass matrix, due to the replica symmetry of Eq. (57), can be put in the form (80) by a proper symmetrization 
of indeces in Eq. (155). The parameters entering in that expression are given by 



3)Ac(l,4)] + iAc(l,4)Ac(2,3) 



Afi(l,2;3,4) = ^M^gli) _ 1 [^(l,4)Ac(2,3) + <S(2, 
.9(1,2) 2p 

Af2(l, 2; 3, 4) = 4) Qj(2, 3) - Ac(2, 3)) - h{2, 3) Qj(l, 4) - Ac(l, 4) 

M3(l,2;3,4) = 15(1,4)5(2,3) 
where we have introduced the notation 

Ac(l,2) = c(l,2)-S(l,2) . 



(156) 



(157) 



Note that we did not symmetrize Eq. (155) over spatial indeces (i.e. over exchanges 1 o 2 and 3 <-> 4). This has not 
to be done explicitly because in any case we are going to apply these operators to symmetric functions only. 
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C. Expression of A in HNC 



Plugging Eq. (153) into Eqs. (97) we get the following two expressions for wi and W2 



1 



W2 



Wi 



1 



fco(l, 2)A:o(3, 4)fco(5, 6)5(1, 3)(5(2, A)S{1, 5)(5(2, 6) 



1,2;3,4;5,6 



1 



PV(1,2) 



(/ * ko)^ Jl,2,3A,5,6 



fco(l, 2)fco(3, 4)fco(5, 6) r(2, 3)r(4, 5)(r(6, 1) - 3f(6, 1)) + (3r(2, 3) - f(2, 3))f(4, 5)f(6, 1) 



(158) 



where we have denoted r(l,2) and f(l,2) respectively the diagonal and off diagonal part of the matrix (151). It 
follows that the exponent parameter is given by 



/d, 



1 fAr^irl 



j^IjM^kliq) [l-pAc{q)Y 



(159) 



D. Computation of ^ and a in HNC 



We will now compute the coefficients of the mass matrix {fi, a, m2, m^) in the HNC approximation. The contraction 

(v) 

of the operator with the zero mode gives 



mi(p) 



dqdk 



(27r 



ko{q)M[P\q,k)ko{k) 



1 

1 



dr 



dq 



MM _ J 

g{r) 2pJ (27r) 



kliq) 



Ac 



Ac 



dq 



(2^) 



D 



fc2(q)Ac 



q] Ac 



so the the fi coefficient defined in Eq. (120) is given by 



fi ~ lim 



dmi {p = 0) 



In order to compute this derivative we recall that from Eq. (74) 

-dg{r) 



nkoir) — lim y/e- 

e-i-o de 



~g{r, e) = ~g{r, 0) + 2^eKkoir) + 0(e) 



Then we can use the replicated Ornstein-Zernike relation [47] to obtain 

Ac(g) 



1 - pAc(g) 



from which 



ciq, e) = c(g, 0) + ^/ico(q) + 0(e) 
Co(<z) = hm ^ = 2Kko{q) [1 - pAc{q)f 

Using these expressions we arrive to the final form for the coefficient p: 



/ dr 



2k 
P 



'^'^ kl[q)[l-pAc{q)f 



(2^) 



In an analogous way we can obtain the expression for the coefficient a 

dmi(p) 



tr = lim 

e->0 dp^ 



p=0 



(160) 



(161) 



(162) 



(163) 



(164) 



(165) 



(166) 
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To compute this expression let us use the foUowing relation 



/ 



P 1 
2+' 



4 g2 



so that the final expression for a is given by 



dq 



1 



8pJ (27r) 



j^kliq) [pAc{q) - 1] 



Ac"{q) 



4 q 
Ac'iq) 



cos^ 61 - 



4 g3 



1 



■fcg(g) (Ac'(g))^os2 



Sy (2^)^ 

where is the angle between the D-dimensional vector q and one of the coordinate axis. In Z3 = 3 we get 

1 ( 1 

= 48^ X '^'i''o(<l)[-iP^<l) - 1] ['?'Ac"(g) + 2gAc'(g)] - g2 {Ac\q)f 

From Eq. (117) we can compute also the parameters m2 and 7713. The expressions for these two quantities are 



(167) 

(168) 
(169) 

(170) 



7712 = 



dq 



(27r) 



^k^{q)£{q) 



- - Ac(9) 
P 



TO3 



dg ,.2,„^,2 
(2^)^ 



(171) 



E. Summary of the numerical results 



Using the HNC approximation for the free-energy, i.e. neglecting 2PI diagrams in Eq. (51), leads to a self-consistent 

(2) 

equation for the order parameter p^^^ via Eq. (56). With the choice of a replica-symmetric structure Eq. (57), and 
setting TO = 1, we obtain a self-consistent equation for the diagonal correlations that coincides with the liquid HNC 
equation, and a self-consistent equation for the off-diagonal parts that takes the diagonal ones as input (see [47] for 
details). Note that although this calculation is possible in any D, here we restrict to D = 3 for simplicity. 

The equations are solved numerically using a standard Picard iteration scheme. We first focus on the diagonal 
part. We start from a very low density p sa 0.2 and gradually increase the density while following the evolution of 
the solution. At high enough density, if we solve the off-diagonal equation starting from a suitable guess for the 
off-diagonal c, we obtain a non-trivial solution. When a non-trivial c has been obtained for a given p, we gradually 
lower the density while following the evolution of c, in order to get very close to the critical point pd where the solution 
disappears. 

For the obtained values of p > pd,^e then numerically compute the derivative of the order parameter g with respect 
to density, both in Fourier and real space. Recalling that this derivative is divergent at the transition, as shown in 
Eq. (74), we determine the precise value of pd by imposing that for fc « 27r in Fourier space and r = in real space, 
the derivative of the correlation function scale as using pd as a fitting parameter. We used these particular 

values of r and fc, since we observed that they were the most sensitive to density changes. The prefactor of this 
square-root behavior is by definition the zero- mode ko in Eq. (74). We check the validity of this scaling by computing 
ko at two different densities slightly above pd, and by checking that fcp does not depend on p close enough to pd- Note 
that here we included the constant k in fep. In fact it is easy to see that the overall normalization of kg does not affect 
any of the physical observables. We show in Figure 1 the typical shape of the zero-mode ko in Fourier and real space, 
the off-diagonal pair correlation function g in real space, and its Fourier transform, all computed at pd- 

We applied this numerical scheme to several benchmark models of monodisperse systems of three-dimensional 
spherical particles that have pair-wise interactions. The models we used are: 

• Hard spheres (HS): 

• Harmonic spheres (HarmS): 

-'•/••»>" (173, 

for various temperatures (note that for /3 — > cxj this potential reduces to the hard-spheres potential), 
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FIG. 1: OfT-diagonal correlations and the zero mode, all computed at pd for hard spheres in D = 3. 



Lennard- Jones (LJ): 

Weeks-Chandler- Andersen (WCA): 

v(r) 

Soft-spheres (SS): 
for n = 6,9, 12. 



v{r) — 4e 



\ r 



4e 



(7) 



12 



-) 

r / 



6 



= ^ (7) 



(174) 



(175) 



(176) 



The lengths are computed in units of the particle diameter rg , and the temperatures in unit of the energy scale of the 
potential e (except for HS, for which temperature is irrelevant). Note that for SS temperature and density are not 
independent variables, so we used the density as a control parameter. 

Once fcoj c and c have been determined at the critical point, we can readily compute the different parameters 
involved in the calculation of the exponent parameter, the prefactor of the correlation length, the prcfactor of the 
Ginzburg criterium, as well as the prefactor of the divergent part of the four-point correlation function. Note however 
that the latter depends on the function / that we choose in Eq. (37) to define our order parameter. In the numerical 
computation we used a box function 



D 



f{x)^{2Ar^\[e{A' 



(177) 
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System 


T 


Pd 


Wl 


■W2 


7712 


m-j, 


(T 




A 


Co 


Go 


Gi 


SS-6 


1 


6.691 


0.121 


0.0845 


-0.229 


0.0273 


0.0484 


0.130 


0.697 


0.601 


224 


0.370 


SS-9 


1 


2.912 


2.41 


1.70 


-1.34 


0.157 


0.405 


1.35 


0.705 


0.548 


34.3 


0.166 


SS-12 


1 


2.057 


8.58 


6.08 


-2.89 


0.328 


0.938 


3.77 


0.709 


0.498 


14.2 


0.154 


LJ 


0.7 


1.407 


33.1 


23.5 


-6.39 


0.719 


2.45 


10.3 


0.709 


0.489 


6.00 


0.108 


HarmS 


10"^ 


1.335 


40.4 


29.1 


-8.34 


0.850 


1.92 


19.3 


0.719 


0.315 


2.82 


0.535 




10"* 


1.196 


51.5 


38.9 


-10.0 


0.957 


2.03 


27.0 


0.756 


0.274 


1.69 


0.622 




10"^ 


1.170 


54.3 


41.5 


-10.3 


0.979 


2.09 


27.1 


0.764 


0.278 


1.66 


0.593 


HS 




1.169 


54.5 


41.5 


-10.3 


0.984 


2.10 


26.7 


0.761 


0.280 


1.67 


0.606 



TABLE I: Numerical values of the coefficients of the effective action and the physical quantities from the HNC approximation. 
For each potential, lengths are given in units of ro and energies in units of e, with ks = 1. Data at fixed temperature, using 
density as a control parameter with € = ps_ — p- 



System p 


Td Wl W2 m2 ms a p. 


A ^0 Go Gi 


LJ 1.2 
LJ 1.27 
LJ 1.4 
WCA 1.2 
WCA 1.4 


0.336 58.2 41.4 -8.94 0.999 3.65 14.2 
1.270 47.9 33.8 -7.97 0.917 3.18 11.1 
0.684 33.7 23.9 -6.46 0.726 2.49 7.24 
0.325 61.0 42.9 -9.65 1.06 3.29 15.1 
0.692 34.6 24.2 -6.69 0.746 2.39 7.21 


0.712 0.507 4.56 0.0937 
0.705 0.536 5.74 0.102 
0.710 0.586 8.52 0.107 
0.703 0.467 4.37 0.180 
0.701 0.576 8.67 0.144 



TABLE IL Numerical values of the coefficients of the effective action and the physical quantities from the HNC approximation. 
For each potential, lengths are given in units of ro and energies in units of e, with fcs = 1. Data at fixed density, using 
temperature as a control parameter with e = Td — T. 



where Q{x) is the Heaviside step function. We choose A = O.Ij'o- For all systems, we give the four main results of this 
paper in Tables VII E and VII E: 

• The value of A given by Eq. (159). 

• The prefactor = ^/^TJ^ of the correlation function ^ = ^qc"^^'', see Eq. (126). 

• The prefactor of the divergent part of the four-point correlation function Gth as given in Eq. (127). 

• The prefactor of the Ginzburg criterion given by Eqs. (148) and (132). 

Moreover let us recall that the quantities reported in the tables are not exactly the quantities defined in the previous 
sections. In particular the coefficients fi and a and the two masses m2 and are computed by including in fcg the 
normalization factor k (hence they are multiplied by k^) and moreover the coefficients wi and W2 are the ones given 
by (97) but again including k in fco (hence they are multiplied by k^) and without the f*ko factor which clearly does 
not enter in the computation of A. In any case the values of the four physical quantities A, Gq, Gi, are independent 
of these normalizations. 

We found that, for all these systems, the dependance of the results on the real and Fourier space cut-offs that are 
needed for the numerical calculation is significant, altering the results on our predictions to within 10"'^. For the case 
of hard-spheres, we checked that the results become increasingly stable when diminishing both cut-offs simultaneously. 
The error that we make because of the finite cutoffs arc however not very important because numerical simulations 
and experiments can not, for the moment, provide results with a better accuracy, due to the difficulty in accessing 
the critical region close to the glass transition. Our results for A are compatible with the ones obtained from Mode- 
Coupling Theory [2, 50, 51], which have also been compared with numerical simulations finding good agreement. 
Unfortunately, not many numerical data for the behavior of the thermal correlation in the /? regime are not available. 
We have checked that our results are roughly consistent with the results of [52] , but more precise simulations would 
be very useful to test our predictions. 

VIII. CONCLUSIONS 

A complete characterization of dynamical heterogeneities in the 13 regime has been obtained using a static equilib- 
rium approach in the framework of the replicated liquid theory. The criticality of the four point density correlation 
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functions has been analyzed through the computation of the stabihty operator of the rephcated Gibbs free energy 
around the glassy solution. This kernel operator has a soft mode that is responsible for the growth of the various type 
of susceptibilities at the dynamical transition. Having identified the soft mode, we have produced a gradient expansion 
for the field theory which describes the fiuctuations of the two point density field that are along the zero mode itself. 
In this way we can focus on the critical part of the quantities we are interested in. Then we have studied the theory 
at the Gaussian level and we have performed a one loop analysis in order to see where the mean field regime breaks. 
This results in a Ginzburg criterion for the dynamical transition. Our approach is completely general and it relies 
only on the fact that the glassy phenomenology can be captured by the one-step replica symmetry breaking scheme. 
Moreover, to produce some quantitative predictions we have computed in the HNC approximation all the observables 
we are interested in: the correlation length, the thermal four points correlation function, the Mode-Coupling exponent 
parameter and all the couplings of the effective replica field theory. 

Our calculations can be straightforwardly extended to other more accurate approximation schemes and this is one 
of the points that must be explored. The question about how to attack systematically the problem in the a regime 
remains open. One might think that replicas are useful mostly in the (3 regime where one is exploring the interior of one 
metastable state, while the study of the barrier-crossing a regime requires a full dynamical approach. However recently 
it has been discussed how replicas can be used to obtain the full long time reparamctrization invariant dynamics [53] . 
This is surprising at first sight because it gives a recipe to obtain some results for the dynamics directly from the 
statics; however the time sector that it explores is the reparamctrization invariant one where quasi-equilibrium holds. 
It would be interesting to see what is the insight that can be gained for dynamical heterogeneities in the a regime 
from the application of this line of reasoning to the replicated liquid case. 

Finally, a recent numerical investigation confirmed the existence of an upper critical dimension du = S for glassy 
dynamics [54]. A more compete comparison between the theoretical predictions obtained here and the numerical 
results would be very helpful to understand the nature of the corrections to mean field that become important 
below d = 8. 
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Appendix A: Double Counting Problem on the ip"^ theory 



Let us consider the (p"^ theory defined by the action in Eq. (10). The generating functional is defined by Eq. (9). 
At the mean field level we have rMF[0] = S[(j)\. Including one loop corrections, we obtain the following expression: 



where K is the following operator 



riL[</'] = rMF[0] + ^Triogif 



Kix, y) = 5{x - y){ml - -[- ^-4,\x)) 



(Al) 



(A2) 



By isolating a free operator Kq(x, y) — S{x — y){m^ — V^) and defining an operator {(f)'^){x, y) — (j)^{x)5{x — y) we can 
rewrite the K operator in the following way 



where ® is the operator (integral) product. Moreover 

Ko\x,y) = 

It follows that 

riL[0] = rMF[</'] + ^Triog 



dp e 



ip{x-y) 



(27r)^ p2 _^ ^2 



1 + 1^^0-'' 



1. 



-Trlog/Co 



(A3) 



(A4) 



(A5) 
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The last term is (/)-independent which means that it can be neglected for our purposes. In fact we want to see what 
happens if we compute the two point function at order g using Til as an action instead of using the bare action 
S. Because we want the two point function at order g we need to perform a Taylor expansion of the extra term 
Trlog[l + {g/2)KQ^ (g) (fp] which is given by 



Tr log 



1 



Tr 



0(ff') = |l?iK) / Axcj?{x) + 0{g^) 



(A6) 



where 



dp 



(A7) 



We now use Til as the bare action, 5'il[v] = TiLiv] ^^'^ compute the two point function. The diagrammatic rules 
are the standard ones for the (p^ theory [25]: 



G(p) 



n. 



(A8) 



where the second diagram is the one originated by the term in Eq. (A6). Note that the second diagram has exactly 
the same expression as the last one. This is the expression of the double counting problem. 

At this point we can compute the Ginzburg Criterion starting by iSiLif/j] = riL[v3]. The action at order g is then 
given by 



1 



dxip{x) —V + ttiq + ^Di{mQ) ip{x) + / dxf {x) . 



(A9) 



It should be clear at this point that the only difference with the bare action S[ip] is a change of the bare mass, 



mo — >■ rriQ + ^Di{mQ). However, we have seen in Section II D that the final expression of the Ginzburg criterion is 



expressed in terms of the renormalized mass only, and is therefore not affected by a change of the bare mass. We 
conclude that whatever microscopic action we use - provided it can be developed in powers of at small tp, which is 
the crucial assumption of mean field theory - will give the same results for the Ginzburg criterion. In this sense the 
Ginzburg criterion can be thought as a check a posteriori of this assumption. 
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